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The Single Server Queue and the Storage Model: 
Large Deviations and Fixed Points 

Moez DRAIEF * 


Abstract 

We consider the coupling of a single server queue and a storage model defined as a 
Queue/Store model in Eld We establish that if the input variables, arrivals at the queue 
and store, satisfy large deviations principles and are linked through an exponential tilting 
then the output variables (departures from each system) satisfy large deviations princi¬ 
ples with the same rate function. This generalizes to the context of large deviations the 
extension of Burke’s Theorem derived in 0C| 

Keywords: Single server queue, storage model, large deviations, Burke’s theorem. 


1 Introduction 

A celebrated theorem of Burke |X asserts that when a Poisson process { a n , n £ Z} with mean 
inter-arrival time 1/A is input to a single server queue whose service times {s n , n £ Z} are 
i.i.d. exponentials (./M/l/oo following the Kendall nomenclature) with mean 1/p, < 1/A, the 
equilibrium departure process {d n , n £ Z} from the queue is also a mean 1/A Poisson process. 
In this sense the mean 1/A Poisson process is a fixed point for the ./M/l/oo operator. In [01E], 
we prove a generalization of Burke’s theorem to the couple of departures and the sequence of 
times spent at the very back of the queue {r n , n £ Z}. The sequence {r n , n £ Z} can also be 
seen as the sequence of departures from a storage model, where s n is the amount of P Supplied 
at slot ra+1, and a n is the amount of P Asked for at the same slot. More precisely, we show that 
under the conditions of Burke’s theorem the couple of output variables from the Queue/Store 
model {(d n ,r n ), n £ Z} has the same law as the couple of input variables {( a n ,s n ), n £ Z}. 

In this paper, we consider the fixed point question at large deviations scaling. A similar result 
has been presented by Ganesh et al. in m for discrete-time queues. Assuming that the input 
variables satisfy a large deviations principle, we show that if their laws are linked through an 
exponential tilting, then the large deviations principle is preserved by the Queue/Store system 
(i.e. the output variables satisfy large deviations principles with the same rate functions as the 
input variables). 

This paper is organized as follows: section [3 is devoted to giving some background on large 
deviations that will be useful for our analysis. The Queue/Store model is presented in section 
01 Section 01 focuses on the workload process. We define the concept of effective bandwidth in 
section 01 Finally, in section 01 we derive a large deviations principle for the output variables 
of the Queue/Store model 1 )1(111) and we give the condition on the rate functions of the input 
variables which are preserved by the Queue/Store operator. 
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2 Large deviations 


Let (Ll,T,P) be a probability space and x = {x n ,n £ N} be a sequence of i.i.d. random 
variables. We define the sequence of partial sums X = {X n ,n £ N} by 

n 

X n — ^ ^ Xi . 

i= 1 

If E|xo| < oo, then X satisfies the strong law of large numbers , i.e., 

lim —- = Exo, a.s. 

ii—xx) n 

We also focus on the fluctuation of a random variable around its mean. Suppose Exq < oo then 
the central limit theorem gives the fluctuations of the scale of 0{\/y/n), of X n /n around Exo- 
More precisely the sequence — Exo^ converges in law a Gaussian random variable with 

mean zero and variance the one of xo- Large deviations theory deals with the fluctuations of 
the scale of 0(1). 

Definition 2.1. Let X be a real Hausdorff space. A function I : X —» M+ = M + U{oo} is a rate 
function if I is lower semi-continuous, i.e. the sets {x : /(x) < a} are closed, for all a £ R. 
In addition, if these sets are compact then I is a good rate function. 


Definition 2.2. A sequence { x n , n £ N} in X satisfies a large deviations principle with the 
rate function I : X —> M +; if for each Borel subset 23 ofX, 

— inf I(x) < lim inf — logP(x n £ 23) 

xG'B 0 n—>oo n 

< lim sup — logP(x n £ 23) < — inf(/(x) , 

n — >oo xG'B 

where 23° denotes the interior of 23 and 23 its closure. 


Throughout the rest of this section, we are interested in real-valued random variables. 

Definition 2.3. The cumulant generating function of a real-valued random variable x is given 
by 

A x (6) = logEe 01 , V0 £ M , 

which can be infinite. 


Using Holder’s inequality and Fatou’s lemma, we check that Ax is convex, lower semi- 
continuous and for all 6 in the effective domain of Ax (i.e. {9 \ A x{9) < oo}), we have 


Aa-W 


E(X 0 e ex °) 

e A x(0) 


We define the Legendre transform of Ax as 


Ix(x) = sup{0x - Ax(0)} • 


(1) 



The function I\ is positive, convex, and lower semi-continuous. If Ax is finite at the origin 
then 

A'y(O) = EXo, and Ix{ EXo) = 0 . 

Moreover, we easily check that 

Ix{x) = sup{0x — Ax(#)} for x > Exo 
6>0 

Ix{x) = sup{6»x — Ax(0)} for x < Exo . (2) 

0<O 

We recall the statement of Cramer’s theorem for real-valued random variables. 

Theorem 2.4. Let {x n , n £ N} be an i.i.d. sequence of real-valued random variables and 
let Ax be its cumulant generating function which we suppose finite at the origin. Then the 
sequence of random variables {X n /n, n £ N*} ; where X n = X 4 L 1 x i> satisfies a large deviations 
principle with rate function lx defined in O): for each closet subset F o/M, 

limsup — log P(S n £ F) < — inf lx (x) . (3) 

n —>oc n x£F 

and for each open subset O of R 

liminf — logP(S n £ O) > — inf Ix(x) . (4) 

n —>oo n xeO 

The inequalities 0 and 0 are known as the upper and lower bounds of the large deviations 
principle. For the different extensions of Cramer theorem we refer to J|. We can have a large 
deviations principle for a couple of independent random variables. 

Theorem 2.5. Let X and y be two real Hausdorff spaces, and { x n , n £ N} and {y n , n £ N} 
be two sequences satisfying large deviations principles on X and y, with rate functions lx and 
Iy■ Suppose that they are independent, then the sequence {(x n ,y n ), n £ N} satisfies a large 
deviations principle on X x y with rate function 

I X ,y(x, y) = I x (x) + I Y (y) . 


In this paper, we will establish large deviations principles relying on indirect methods. Once we 
have a large deviations principle for one sequence of random variables, we can effortlessly obtain 
large deviations principles for a whole class of random sequences, namely those obtained via 
continuous transformations. We present the contraction principle, the tool that enables this. 

Theorem 2.6. If {x n , n £ N} satisfies a large deviations principle on X with rate function lx, 
and if f : X —> y is a continuous function, then the sequence {y n , n £ N} defined by y n = f(x n ) 
satisfies a large deviations principle on y with rate function 

Iy(y ) = inf Ix(x) . 

X, f(x)=y 

Sometimes / is “almost continuous”, i.e. the sequence {y n , n £ N} is close to a sequence 
{/(x„), n £ N}, where / is continuous. 
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Definition-Proposition 2.7. Let {x n , n £ N} and {y n , n £ N} be two sequences on X. They 
are said exponentially equivalent, if for all e > 0, 

limsuplogP(||x n - y n \\ > e) = -oo . 

n—xx) 

In this case, if {x n , n £ N} satisfies a large deviation principle then {y n , n £ N} satisfies a 
large deviations principle, and lx = Iy ■ 


Theorem 2.8. If {x n , n £ N} satisfies a large deviations principle on X with rate function lx, 
and if {y n , n £ N} and {f(x n ), n £ N} are exponentially equivalent, with f : X —> ^ continuous, 
then {y n , n £ N} satisfies a large deviations principle on ^ with rate function 


Iy{y) 


inf I x {x) . 

xeX, f(x)=y 


For a thorough presentation of large deviations’ results and applications, we refer to p]. 


3 The Model 


Let A = {A n , re £ Z} be a point process and assume that Aq < 0 < A\ and A n < A n+ 1 , Vn £ 
Z. We define the Mi'j.-valued sequence of r.v’s a = { a n , re £ Z} by a n = A n+ \ — A n . Let 
s = {s n , n £ Z} be another Ml-valued sequence of r.v’s. The sequences a and s are the input 
variables of the model. 


Define the sequence of r.v.’s T> = {D n , n £ Z} by 


D 


n 


n 


sup 

k<n 


Ak + ^2 Si 

i=k 


(5) 


A priori the D n ' s are valued in ]R U {+oo}. Assume that a and s are such that the D ra ’s are 
almost surely finite. Set d n = D n+ \ — D n and d = {d n , n £ Z}. Define an additional sequence 
of r.v.’s r = {r n , n £ Z}, valued in M?j_, by 

r n = min(D n , A n+1 ) - A n . (6) 


The sequences d and r are the output variables of the model. In view of the future analysis, it 
is convenient to introduce the sequence w = {w n , n £ Z} of random variables valued in M + , 
defined by 


n— 1 


Wn = D n - s n - A n = sup 

k<n— 


J£' 


Si CLi 


These random variables satisfy the following recursion (Lindley’s equation) 


(7) 


w n +1 — \w n H - s n a. n ] 


( 8 ) 


Using the variables w n , we can give alternative definitions of D n and r n \ 


Ml < n, D n = 


wi 


+ A-i+ . 


i=l 


V max 

l<k<n 


Ak +< 


(9) 


i=k 


r n = min+ s n , a n } = s n + w n - w n+ 1 . (10) 

We now interpret the variables defined above in two different contexts: a queueing model and 
a storage model. 


4 









3.1 The single-server queue 

We are concerned with a single server queue where each customer is characterized by an instant 
of arrival in the queue and a service demand. Customers are served upon their arrival in the 
queue and in their order of arrival. Since there is a single server, a customer may have to wait 
in a buffer before the beginning of its service. Using Kendall’s nomenclature, our model is a 
././1/oo/FIFO queue. The customers are numbered by Z according to their order of arrival 
in the queue (customer 1 being the first one to arrive strictly after instant 0). Let A n be the 
instant of Arrival of customer n and s n its Service time. Then the variables defined in ©-© 
have the following interpretations: 

• D n is the instant of departure of customer n from the queue, after completion of its 
service; { d n , n G Z} is the sequence of inter-departure times; 

• w n is the waiting time of customer n in the buffer between its arrival and the beginning 
of its service; 

• r n is the time spent by customer n at the very back of the queue. 

The variables {r n , n G Z} are less classical in queueing theory ;Q:. 

3.2 The storage model 

Some product P is supplied, sold and stocked in a store in the following way: Events occur at 
integer-valued epochs, called slots. At each slot, an amount of P is supplied and an amount of 
P is asked for by potential buyers. The rule is to meet all demands, if possible. The demand 
of a given slot which is not met is lost. The supply of a given slot which is not sold is not lost 
and is stocked for future consideration. 

Let s n be the amount of P Supplied at slot n + 1, and let a n be the amount of P Asked for at 
the same slot. In this context, the variables in ©-© can be interpreted as follows: 

• w n is the level of the stock at the end of slot n. It evolves according to ®; 

• r n is the demand met at slot n + 1, see equation GDI; it is the amount of P departing at 
slot n + 1; 

The variables {D n , n G Z} do not have a natural interpretation in this model. 

It is important to remark that while the equations driving the single server queue and the 
storage model are exactly the same, the relevant variables are different. The important variables 
are the ones corresponding to the departures from the system. The departures are coded in 
the variables {d n , n G Z} for the single server queue and in the variables {r n , n G Z} for the 
storage model. For a more detailed discussion of these two models we refer to j3|. 

3.3 Rare events 

We recall that w n is the waiting time of customer n before it starts its service (respectively the 
level of the stock at the end of slot n) and it is given by 


n— 1 



(ii) 
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We assume that the sequences a = { a n , n E Z} and s = {s n , n E Z} are stationary and 
ergodic. Under the stability condition Eso < Eao, we can use the above expression of w n 
to give its asymptotic behaviour using the large deviations properties of the input variables 
a = { a n , n E Z} and s = {s n , n E Z}. We start with an example. 


Example Suppose that customers arrive in a deterministic fashion, i.e. a n = 1, Vn, and that 
the sequence {s n , n E Z} is i.i.d. with 

P(s n = 2) = 1 - P(s n = 0) = p < i . 

The process {run, n E Z} is a discrete time birth-and-death process with stationary distribution 

p <»o > i)=(j-y ■ 

We get 

where 5 = log 

This approximation remains valid under general conditions on the input variables, with an 
expression of S that depends upon the rate functions associated to the sequences a and s. 
Indeed it has been proved El El CS! that the stationary version of w n satisfies the following 
property 

lim — logP(u;o > nq) = — Sq , (12) 

n—>oo 77, 

or alternatively 

¥(wo/n > q) x e~ nSq . 

First we give a heuristic proof of this result which will be asserted rigorously in section 0J We 
assume the sequences {a n , n E Z} and {s n , n E Z} i.i.d. mutually independent with 

A a {6) = logEe eao , A s = logEe es ° , 

finite near the origin. Let x n = s n — a n for all n E Z, and Xk = ^i=i x -i for all k E N then 
wo = sup k>0 X k , with Xo = 0. Using Cramer’s theorem (Theorem 12.411 . the sequences {a n , n E 
Z} and {s n , n E Z} satisfy large deviations principles on M with rate function I a and Is the 
Legendre transforms of A^ and As. Thus the sequence {x n , n E Z} satisfies a large deviations 
principle with rate function lx, which is the Legendre transform of A x(Q) = A g(6) + Aa(—9), 
i.e. P(^- > x) x e ~ nI x(x) Moreover, 

P(w 0 > nq) = P(supX fc > nq) = P(U k > 0 {X k > qn}) 
k> o 

OO 

< 5^P (X k >q n ). 

k =1 

Since P(X^ > nq) = P(^r > x e~ kIx ^ nq ^ k \ we check that 

P(u;o > nq) x e~ nq n ^l k . (13) 

k=\ 


— logP(u;o > nq) = —5q , 
n 
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We conclude using the principle of the largest term, i.e. 


UA iyO.g/fc) .. 

\ g ^ nq/k g 

k=1 

where 5 = inf x> o x x and we get ea>- The principle of the largest term translates the fact 
that rare events occur in the most likely way. Indeed, the dominant term in <tm which gives 
the explosion of the waiting times in the queueing system (or the overflow of the stock in the 
storage model), happens following the most probable scenario. 


4 Large deviations for the workload process 

First, we give the assumptions under which the approximation m is fulfilled. For more general 
assumptions we refer to the paper by Ganesh et al. [El- 


Assumptions: 

(i) The sequences {a n , n G Z} and {s n , n £ Z} are i.i.d., mutually independent. Their 
cumulant generating functions are given by 

Aa{0) = log , As(<9) = log Ee 0s ° 

We assume that both and A 5 are differentiable near the origin. 

(ii) The stability condition 


l/7i = A^(0) = Esq < Eao = A^(0) = 1/A 


is satisfied. 


Under ( i ) and (ii), we have 

Proposition 4.1. The sequence {wo/n, n G N*} satisfies a large deviations principle with good 
rate function I\v(q) = <*><?, he. linr, woo ^log¥(wo/n > q) = — 5q, with 

6 = inf A - s (- i__ _ gU p . \ s (0) 4 - A a(—Q) A 0} . (14) 

0 <a<s s — a 


Proof. First we prove <©• Recall that x n = s n — a n , X n = x -i> an d 

w 0 = sup X k . 
k> 0 

By the independence assumption, we have 

A x(0) = A s (0) + Aa(-O) . 

By direct application of the contraction principle, we check that 

I x (x) = inf {I s (y) + Ia(v - z)} ■ 

y>x 
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Let 0 < inf x >o I x (x)/x then 


0 < inf Ix(x)/x 

x>0 


44 

0 < Ix{x)/; 

x,Vx > 0 

44 

0x — I x {x) 

IA 

o 

< 

IV 

o 

44 

sup{ 0 x — Ix(x)} < 0 


#>0 


44 

A x (0) < 0 . 



The last equivalence is due to the fact that Exo < 0 (stability condition) and the equation ©. 
We proved that 

inf ^ A ’ S ( a ' = inf I x (x)/x = sup {0 : A s{0) + AM— 0 ) < 0 } . 

0<a<s S — a x>0 


Lower bound: For q > 0, we have P(u>o > q) > P(X\ > q). Notice that for p > we have 


P (™o >q)> > p) . 


Since - > 1 1 


> - r 7 -i . we get 
q — p \q/p 1 \ ’ B 

liminf - logP(tt> 0 > q) > ~ liminf > p) 

9 ->°o n v >oo \a v\ V o/o / 


p q- 
1 


p n—»oo 77 , 


= - liminf — logP( —X n > p) > — lx(p) , 


n 


for all p > 0. We conclude that 


1 


liminf - logP(u;o > q) > — inf —I x (x) = — inf — inf {Is(y) + Ia(u ~ *)} 

9^00 q x >0 x x>0 x y>x 


= — inf 


Is(s) + Ia((l) 


0<a<s s — a 


• Upper bound: Recall that A x (0) = A g(9) + A a(—9)- By the stability condition, we have 


A^(0) 


log 


E(si) 

E(ai) 


< 0 . 


Moreover, Ax is differentiable near 0 with Ax(0) = 0. Thus, there is a constant 0 > 0 such 
that Ax(0) < 0 and Ee eXn < oo. Applying Chernoff’s bound, we get, for n € N 

P(X n > q) < e~ eq Ee eXn . 


This leads to 


lim sup - logP(X n > q) < —0 . 

q—>oo Q 


For N € N, we check that 


pf max X n > q] < N max P(X n > q) . 
\0<n<N ~ V ~ 0<n<N V ~ 


Allowing q go to infinity, 


lim sup —log pf max X„ > q\ < max lim sup log — P(X„ > q) < —© . 
g—>oo q \0<n<N J 0 <n<N ^ q K 
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(15) 


We need now to have a bound for large values of n. Applying the union bound, we have 


1 sup 

n>N 


x n >q)<Yl >4) < e" 09 Ee 


ex n 


n>N 


n>N 


Since y logEe eX/1 = Ax(0) < 0, there is 0 < e < — Ax(0) et Nq £ N such that 


1 


Vn > Nq, logEe 0Xn < —e . 


n 


Combining m and ©, we obtain 


sup 

u>Nq 


X n >q)< e~ eq J2 ^ < l 


09 

_g — (V©+l)e 


n>N e 


— e 


To sum up, we derived the following bound 

lim sup - 

q^oo q 

Therefore, for 0 > 0 such that Ax(0) < 0, we have 


sup X n > ) < —0 . 

u>Nq 


lim sup - log P ( sup X n > q ) < — 0 . 


q^oo q 


n> 0 


We conclude using m- 


(16) 


□ 


This result can be interpreted naively in terms of the following approximation 

P(sup X n >q)& sup P(X n > q) . 

n n 

This is nothing more than a consequence of the fact that rare events occur in the most likely 
way. Moreover, using this approximation one can predict the frequency at which the waiting 
time (or the level of stock) overflows a given threshold. Indeed if (II21) is fulfilled, then we 
represent the graph (in logarithmic squale) of the frequency of exceeding a level q. The above 
approximation (Proposition 14.11) inspired numerous applications (described by Courcoubetis et 
al. [2) of large deviations to the analysis of the statistics of networks of queues. 


5 Effective bandwidth 


The bandwidth of a traffic ffow has become a classical feature in communication networks 
literature. We are interested in the notion of effective bandwidth introduced by Kelly in ESI, 
which gives an analytical way of describing the properties of a stochastic flow by means of rare 
events occurring in the network through which it passes. 


In practice, the queues and stores we are interested in have finite capacities, i.e. the queue 
(respectively the store) rejects customers (respectively stock) each time the waiting time (re¬ 
spectively the amount of stock) overpasses a given threshold. In the queueing context, we 
assume {s n , n £ Z} given and we seek the (deterministic) minimal value a p of inter-arrival 
times such that the probability that the waiting times exceeds the threshold is less than a fixed 
value p. For the storage model, we assume { a n , n £ Z} given and we seek the maximal value 
s p of product arriving at the store at each slot of time, such that the store rejects it with a 
probability less then p. More precisely, we want to identify the arrivals to both models such 
that 

P(u>o>q)<P, (17) 


for given values of q and p. 
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Queue Let a n = a, Vn E Z, i.e. Aa(9) = aO. If s is i.i.d., then by Proposition 14.11 

P(u >0 > q) < e- 5{a)q 

with 

5(a) = sup{0 : 9 > 0 and A g(9) < a9} . 

The minimal inter-arrival time for the inequality m to be fulfilled is given by a p such that 

a p = inf{a : a > 0 et e~ s ^ q < p} . 

Let 9 P = — then 5(a p ) = 9 P and 

A s(8p) 

a P = — ft -• 

Up 

The variable a v is the effective bandwidth of the queue. For more details we refer to Kelly’s 
review Ill- 


Store Let s n = s, Vn E Z, then A s(9) = 9s and 

5{s) = — inf{0 < 0 : A a(9) < s#} . 

the inequality m is satisfied if no more than s p amount of stock arrives to the store at each 
slot with 

s p = sup{s > 0 : e~ 5 ^ q < p} . 

The quantity s p is the effective bandwidth of the store, with 

Aa{ 9 p ) 


We generalize this notion to non-deterministic arrivals by introducing the functions ola{9) and 
a s (9) given by 


oia(8) = 


A s(8) 


as(9) = 


A A (0) 


9 1 9 ’ 

representing the effective bandwidths of the queue and the store, respectively. 


6 Principal of large deviations for the output variables 

We proved, in section 01 a large deviations principle for the workload process using classical 
techniques. It is generally hard to apply these techniques to derive large deviations principles 
for the output variables. In this section, we will apply the contraction principle to this end. 
Following the outline of the proof of the existence of fixed points for discrete-time queues in 
El, we use theoretical results on large deviations for continuous-time processes. 

First, we go back to the workload process to illustrate this method. For n E N*, we define 

n n 

-A-n — ^ ^ tt-ii S n — ^ ^ S—i . 

2=1 2=1 
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It is obvious that wq = swp k>0 (S k — A k ). 

To apply the contraction principle it is crucial to define an adapted topology under which the 
workload process is obtained through a continuous mapping of the input variables. Therefore, 
we define the polygonal approximation in n of a given sequence {X k , k G N*} by 

X n (t) = —X^ nt j + (t — ~~)(^L ri tl+ 1 — t> 0 . 

Let A n and S n the polygonal approximations in n of the sequences {A k , k G N} and {S k , k G N}. 
Notice that for rGN’, 

— = - sup (S k - A k ) = - sup(5 LntJ - Ay nt j) 
n n km n t > 0 L J L J 

= sup (S n (t) - A n (t)) . 
t> o 


A function x is absolutely continuous (on M) if for every e > 0, there is rj > 0 such that for 
any finitely many disjoint open intervals (a*, bf), i = 1. .. n, satisfying Y^i=\(Pi ~ &i) < we 
have X)iLi(/(&i) — f( a i )) < e . If x is absolutely continuous, then its derivative x' exists almost 
everywhere and we can write 

p V 

x(v) — x(u) = / x(t)dt. 

J U 


For fi > 0, we define (respectively the set of continuous functions (respectively absolutely 
continuous) x : M + —> M with x( 0 ) = 0 and 


lim LL 

t — > + 0O t “|“ 1 


l/fl < oo , 


equipped with the norm 


sup 

teR + 


x(t) 

t + 1 


We now focus on large deviations on processes in continuous time in 


(18) 


Definition 6.1. A sequence of processes {X n , n G N}, where X n G satisfies a functional 
large deviations principle with linear geodesics, with instantaneoiLS rate function I, if 

(?') the function I is a rate function, with 1(1/n) = 0 , 

(ii) the sequence {X n , n G N} satisfies a large deviations principle with rate function 

3 x ((j)) = I ^ +OC W) dt */ ^ G A u 

1 +oo otherwise. 


Dernbo and Zajic explored extensions 0| of large deviations principles to processes in contin¬ 
uous time under the topology of uniform convergence on compact intervals. However, this 
topology is not appropriate in the context of queueing systems. We introduce a coarser topol¬ 
ogy corresponding to the norm m The following large deviations principle is due to Ganesh 
and O’Connell [IX|. Notice that if x = { x n , n G N} an i.i.d. sequence with mean l//i, 
then for n G N*, X n G C^, where X n is the polygonal approximation in n of the sequence 
X = {X n , n G N} with X n = Ya=i x i 
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Theorem 6.2. m Let x = {x n , n G N} be a sequence of real-valued random variables and 
kx its cumulant generating function, which we suppose differentiable near the origin. Then the 
sequence of polygonal approximations {X n , n G N*} of X = {X n , n G N}, where X n = Y/h=i x i> 
satisfies a functional large deviations principle with linear geodesics on equipped with the 
norm with mean 1/p = A x (0) and instantaneous rate function lx, i.c. 

J x (0) = Jio + °° I x{fi'{t))dt if fie A^ 


+oo 


otherwise. 


Under (i), the processes A n and S n satisfy functional large deviations principles with linear 
geodesics and instantaneous rate functions 4 a and ?s defined by 


rwp iifieA x 

| Too otherwise , 


with I A (x) = sup 0gR {0x — Aa(#)}, the Legendre transform of k A . We define Is and 4s dr the 
same way. 


We check (see M) that / : e M x e A 


defined by 


fifi, fi) = sup - fi{t)) , 

t> o 

is continuous on x G x for the topology induced by the norm (H. Indeed / is not continuous 
under the topology of uniform convergence ( m Example 5.2]). 

Applying the contraction principle, we check that {w\/n, n G Z} satisfies a large deviations 
principle on M + with rate function 

roc 

Iw(q) = inf{ / I At s(fi'(s)) ds I {fii,fi 2 ) eA x x A^, f(fi u fi 2 ) = q} 

Jo 

roo 

= inf{ / I A ,s(fi'(s)) ds | (fii,fi 2 ) G A x x A^, sup (fi 2 (t) - fifft)) = g} , 

Jo t> 0 

where I A> s(fi'(s )) = I A {fi\ (*)) T Is{fi 2 (s)). Combining this with Proposition l4.il we have 
Corollary 6.3. The solution of the following optimization problem 

{ Minimizing Jq 00 I A: s(fi'(s )) ds pour fi = {fi\, fi 2 ) G A x x A^ 
subject to sup t> o{<fa{t) - fiiifi)) = q 

is given by I wifi) = dq where 6 is defined in m 

In the following paragraph, we make use of the same techniques to give large deviations prin¬ 
ciples for the output variables 


dn — On + UT+1 VJ n T Sn+1 S n , r n — S n W n -f-i T w n . 
We define D = { D n , n G N} and R = {R n , n G N} by 

n n 

D n — ^ ^ d—f , R n — ^ ^ r _2 . 


(19) 


i= 1 


i =1 


Let D n and R n be their polygonal approximations in n. Under (i), (ii) and by Loynes’ scheme 
m , we check that 

{D n ,R n ) G Ca X 6(i ■ 
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6.1 Optimization problem 

To make the proof as clear as possible, we introduce a new variable 

bn — Q>n T ^n+1 i 

which stands, in a queueing context, for the time between instances of beginning of services for 
customers n and n + 1 . Let B = {B n , n E N} defined by B n = Y^=i b-i and B n its polygonal 
approximation in n then B n E Under ( i ), we have 

Lemma 6.4. for 0 < t < 1 fixed, the sequences { D n (t ) , n E N*} and {B n {t) , n E N*} are 
exponentially equivalent. 


Proof. Let 7 > 0 , 

P( sup | D n (t) - B n (t) | > 7 ) < 

0<«<1 

< 

< 

Since Implex, log((ro + l)e _7 *' n ) = —00 and thanks to assumption (i) (As ( 7 ) < 00 for 7 close 
to 0 ), we are done. □ 

Thus if one of the sequences {D n (t) , n E N*} and {B n (t) , n E N*} satisfies a large deviations 
principle then the other does too and they both have the same rate function. Since the expres¬ 
sion of B n (in terms of the input variables) is easier than the one of D n , we concentrate on B n 
and R n . In fact, 

Lemma 6.5. The sequences B n (t ) and R n (t ) satisfy 

Bnfit) = A n (t) + swp{S n (s) - A n (s)} - sup {(S n (s) - S n (t )) - (A n (s) - A n (t))} 

s>0 s>t 

Rnfit) = s n (t) - sup{5 n (s) - A n (s)} + sup{(5 n (s) - S n (t)) - (A n (s) - A n (t))} 

s>0 s>t 


P( sup Is_ fc - s 0 | > ray) 

0<fc<n 

n 

£p(i*- fc ~ s °i > 

k =0 

n 

e - 7I /n J- Ee7 |,_ fc -, 0 | < ( n + 1 ) e -7^ e 2 As(7) . 
k =0 


Proof. We check the result for B n (t ), the proof goes along the same arguments for R n (t). Since 
B n (t) is a linear interpolation of B\ nt \/n, we will prove the result for B\ nt \. First, we check 
that 

B\nt\ -d[ni| T Wq W_ \nt\ • 

However, we have u>o = sup s>0 {5' n (s) — A n (s)} then 

w-i nt j = sup{ (S n (s) - S n (t)) - (H„(s) - A n (t))} . 

S>t 


□ 
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Theorem 6.6. Under assumptions ( i ) and (ii), the sequence {{^r, -yy, ^r), n E N*} satisfies 
a large deviations principle on Ml with rate function 

J(xi,x 2 ,w) = inf {5(w - xi + x 2 ) + Ia,s{x 2 , ah); inf g(q, r, vi, v 2 )} , (20) 

o 


where 


g(q,r,v 1 ,v 2 ) = t! a ,s 


x 2 -V 2 X2-V2+W 


+ (1 - t)Ia,h 


V\ v 2 - q 


+ 5q 


-1 — r 1 — r 

with 6 = inf 0<a <s lA< ' a ^ I a s ^ and C = {(q, r, v\, v 2 ) \ q > 0, 0 < r < 1, x 2 -v 2 +w = x\—v\+q}. 


Proof. By Loynes US! and lemma 16.51 we check that (B n , R n ,^) =<f>{A n , S n ) with 

$ : e A x e M ^ e A x x R+. 

Let <f> = ($ 1 , <h 2 , d> 3 ) and 0 = (0i, 0 2 ) € C A x we have 

$1 (0) (*) = Mt) + sup{0 2 (s) - 01 (s)} - sup{(0 2 (s) - 02 (i)) - (01 (s) - 01 (i))} 

s>0 s>t 

$2 ( 0 )(t) = Mt) - SUp( 0 2 (s) - 01 (s)} + SUp{( 0 2 (s) - 0 2 (t)) - (01 (s) - 01 (*))} 

s>0 s>t 

$3( 0 )(*) = SUp{0 2 (s) - 01 (s)} . 

s>0 


( 21 ) 


Applying the contraction principle, the sequence {(B n , R n ,wo/n), n E N*} satisfies a large 
deviations principle with rate function 

roo 

J(-01,02,w) = inf{ / I A ,s(<l>i(s),(f> 2 (s))ds I <f>(01,0 2 ) = (01,0 2 ,w)} . 

■Jo 

This variational problem is generally hard to solve, we restrict ourselves to t = 1. By the law of 
large numbers, {(B n ( 1), R n ( 1)), n E N*} = {(^fq -^f-), n E N*} converges to the mean values of 
the output variables, whereas the large deviations principle gives the fluctuations around these 
values. The sequence {(^fq ^ ,wo/n ), n E N*} satisfies a large deviations principle with rate 
function 

/•OO 

J(x\,X 2 ,w) = inf{ / I J 4,5(0 / i( s )>02( s )) ds I [^(01,02)](1) = (xi,x 2 ,u;)} . (22) 

■Jo 

We introduce the variables <70 and (q given by 

q 0 = sup{0 2 (s) - 0i(s)} 

s>0 

<?i = sup{( 0 2 (s) 0 2 ( 1 )) (01 (s) 01 ( 1 ))} ■ 

S>1 

We rewrite m as follows 

xr = 0 i(l) + q 0 - qi, x 2 = 0 2 ( 1 ) - qo + qi, w = q 0 ■ 

To solve this variational problem we condition on the value of q\ to be equal to q. By the 
stability condition A < /a, we check that q is finite. Then J{x \,x 2 ,w) is the solution of the 
following optimization problem 

rl roo 

j I A ,s(^'{s))ds + J I A ,s(4>'{s))ds over 0 E A\ x q> 0 
qi = q, Xi = 01 (1) + qo - qi, x 2 = 02 ( 1 ) - qo + qi, W = q 0 , 


Minimize 
subject to 
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where 0 = (0i(s), 020 O) and /a,s(0'(s)) = U4(0i( s )) + -Ts( 02 ( s ))- 

Using the auxiliary variable y allows the decomposition of the initial optimization problem into 
two minimization problems on s E [0, 1] and s > 1 coupled through q. More precisely, we 
suppose the variables x±, X 2 and w = go given, if we fix q = q\, then 0i(l) and 02 ( 1 ) are fixed. 
For s > 1, the only constraint on 0i and 02 is qi = q- Consequently, the first minimization 
consists of 

/ OO 

/a,s(0 1 ( 5 ), 02 (s))ds, subject to gi = g . 

The second minimization consists of 

Minimizing / /a,s(0'i( s )> 02( s )) rfs 

Jo 

subject to £1 = 0i(i) + g 0 - <?i, = 02 ( 1 ) - go+ gi, w = qo, 

knowing gi = g. In both cases the minimization is for 0 E A a x A /( . 

For the first problem, let 0 E x A M with 0$(s — 1) = 0j(s) — 0*(1), s > 1, i = 1, 2. By lemma 
16.hi this optimization gives the rate function associated with w$ given in Proposition 14.11 and 
the minimum is 5q where 5 = inf a<s . 

Hence, to determine J(x\,x 2 ,w) it remains to 

Minimize [ lA{<t>i(s)) + Is{<t> 2 (s))ds + 6q , (23) 

00 

for (0i, 02) E A\ x subject to 

' <1 >0, w = sup s > o {0 2 (s) - 0i(s)} 

< Xi = 0i(l) + w - q 
„ x 2 = (j) 2 (l) - w + q. 

• Suppose sup s > o {0 2 (s) -0i(s)} = sup s >i{0 2 (s) -0i(s)} , then w = g +0 2 (1) - 0 2 (1), he. 

£1 = 02 ( 1 ), £ 2 = 0i(1), q = w-x 1 + x 2 . 

The optimisation problem da is fulfilled by linear functions 0i and 0 2 and we have 

5{w -xi + x 2 ) + Ia(x 2 ) + I s (x 1 ). 


If 02(s) — 0i (s) reaches its maximum for r E [0,1] then 

xi = -q + 0i(l) - 0i(r) + 0 2 (r) 

x 2 = q + 0 2 (1) - 02 (t) + 01 ( t ) • 

We define y; = ^ j _ 12 and % E A,\ x A 


70 


Xi(*) = 


Vit if t E]0, r] 

UiT + Zi(t -t) if t E]r, 1] 


Using the convexity of I A and I 5 and Jensen’s inequality, we have 


lA,s(x'i(s),X2( s )) ds < / ^4,s(0i( s )>02( s ))ds 
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Since I A ,s(x'i(s), X 2 ( s )) ds = rI A,s(y) + (1 - r)I A:S {z ) where y = (yi,y 2 ) and z = 
(zi,z 2 ), it remains to minimize 

5q + Tl A>s (y) + (1 - t)I Ai s(z ) , 


subject to 

q > 0, 2/i > 2/2 
xi = -q + (l- t)zi + Ty 2 
x 2 = q + (l — t)z 2 < (1 — t)z\ 
q + (1 - r)z 2 < (1 - t)zi. 

Let v\ = (1 — t)z\ and v 2 = (1 — t)z 2 + q, we check that y\ = X2 ~ V2 ; y 2 = w + x 2 -v 2 anc q 
J(xi,x 2 ,w) = inf{<5(«; - x 1 + x 2 ) + I A ,s(x 2 , xi), inf g(q, r, Vi, v 2 )} , 
g et C defined in Theorem Ki.Cil 


□ 

Let us give an interpretation of this result in a queueing context. Indeed, the most likely ways 
to get departures with mean x\, times spent at the back of queue with mean x 2 and waiting 
times with mean w are when 

• all customers belong to the same busy period (the system is never empty). Thus the 
mean of the departures is equal to the mean of the services and the average time spent 
at the very back of the queue is equal to the mean of the arrivals, 

• or let (q,r,vi,v 2 ) G C the variables where g reaches its minimum. Then the customer 
—nr finds the system empty, and the optimal path splits into two periods. Customer —n 
waits nq time before starting its service. During the first period the average of arrivals 
and services are X2 ~ V2 and X2 ~ V2+W anc j it consists of nr customers. During the second 
period the mean of services is an d the arrivals occur with mean . 

The equation PU|) determines the rate functions of the output variables in terms of the rate 
functions of the input variables. A natural question is when are these functions equal? More 
precisely, when do we have ( Id,Ir ) = (La, Cs)? We start with an example which translates the 
result stated in [?j in terms of large deviations. 


Example: the M/M/1 queue Suppose that {a n , n£Z) and {s n , n G Z} are i.i.d. mutually 
independent and for t > 0, 

P(ao > t) = e~ xt , P(so > t) = . 

For 6 < A < /r, 

AA«) = i°s(Lr«) ■ A s(») = lo «(^) ■ 

Since 6 = sup{0 > 0, A s{0) + A A (—0) < 0}, we check that 5 = g — A. Moreover we have 
I A {a) = A a- log(Aa) - 1 , I s (s) = gs - log (gs) - 1 . 

Resolving the above optimization problem, we have 

J(xi,x 2 ,w) = Aaq - log(Axi) - 1 + gx 2 - \og(gx 2 ) - 1 + (g - X)w 

= I A>s (xi,x 2 ) + 6w . 
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6.2 Existence of fixed points 


Suppose {a n , n G Z} and {s n , n G Z} i.i.d., mutually independent with means 1/A and l/fi. 
They satisfy large deviations principles with rate functions I A and Is differentiable on R + . Let 
0 > 0 such that 

I's(x) - l' A {x) = 0 , Vx G M+ . (24) 

We can rewrite this relation in the following which will be useful for future analysis 

I A (x ) - Ia{v) = Is(x) - Is(y ) - P(x - y ), V(x, y) G (M+) 2 . (25) 

Before stating the fixed point result for the Queue/Store operator for the rate functions, we 
give an interpretation of the relations PH) and (1251) in terms of exponential tilting of measures. 

Definition 6.7. Let 0 > 0, we say that a measure u' is the 0-exponentially-tilted measure of a 
given measure u if for all i6R, 


du' e@ x 

du X fed x dv{x) 


(26) 


where [dv'/dn)(x) denotes the density of v' with respect to u. 


Proposition 6.8. Suppose {a n , n G Z} and {s n , n G Z} satisfy relation 1^41 ) or H3) , with 
marginals u A and vs- Then v A is the 0-exponentially-tilted function of ns- 


Proof. Recall that I A ( 1/A) = 0. We check that 


A a ( 0) = sup[0x — I A (x)} = sup[0x - Is(x) + Is(l/X) + 0{x — 1/A)] 

icGM icER 

= sup[(0 + 0)x — Is(x)] — — I si)-/ A)] 

x£S. A 

= A s (0 + 0) - A s {0) ■ 

Considering the exponentials on both sides, we have 

fm+ e {e+P]x dv s (x) 


I 


f Rt e^ s ( X ) 


(27) 


□ 


Lemma 6.9. Assume that {a n , n G Z} and {s n , n G Z} satisfy assumptions ( i ) and (ii). Then 
the tilt coefficient 0 is given by 

0 = sup{6» : A s (0) + A a (-0 ) < 0} . 
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Proof. Let 6 = — (3 in (1271) . we have 

A a (-/3) = As(0) - A S (J3). 

Since Ag(0) = 0, we obtain A A (—/?) + A g(/3) = 0. Thanks to (1T21) . we have 

/3 < sup{0 > 0, A s{9) + A a(~9) < 0} = S . 

The function G(9) = As(9) + A A (—0) is convex with G(0) = 0. Moreover, by assumption ( ii ), 

G'(0) = A' 5 (0)-A / a (0)<0. 

This means that G is equal to zero twice on [0, /?], at 0 and (3, then 

(3 = sup{0 > 0, A s(9) + A a{~9) < 0} = 6 . 

□ 


Theorem 6.10. Let 6 > 0 and the sequences {a n , n G Z} and {s n , n £ Z} the input variables 
of the model with va the 5-tilted measure of v$ (or ug the (— 5)-tilted measure of u a). Then 
the sequence ^r,wo/n), n € N*} satisfies a large deviations principle on with rate 

function 

J(xi,x 2 ,w) = 6w + I A ,s(x 1 ,x 2 ) ■ (28) 


Moreover, we have 


Id = I A and Ir = Is . 


Proof. By Theorem lfi.61 the sequence {(D n (l), R n (l), w\/n), n £ N*} satisfies a large deviations 
principle on Ri. with rate function 

J(xi,x 2 ,w) = inf{5(u; - aq + x 2 ) + I A ,s(x 2 , aq); inf g(q, r, v±, v 2 )} , 

(J 

We start with the first term of this optimization problem, i.e. 

aq = 02 ( 1 ), x 2 = 0 2 (l) q = w-x i+x 2 . 

We have w = aq — x 2 + q, and thanks to (ESJ, 

6(w - x 1 +x 2 ) + I A ,s(x 2 ,x 1 ) = 6(w - x 1 +x 2 ) + I A ,s(xi,x 2 ) + 6(xi - x 2 ) 

= 6w + I a ,s(x i,x 2 ). 

For the second term, let (q,r,vi,v 2 ) £ C, 


• If V\ = V‘2 , from the convexity property, then 


g(q,r, vi,v 2 ) 


> 


5q + Tl A ,s(y 


x 2 - v 2 x 2 - v 2 + w 


) + (1 - t)I a ,s( 


v 2 


1 - T 


hq + Ia,s{x 2 ,x 2 +w -q) = 5w + I A ,s{Xi,x 2 ) . 


V2 - q \ 

1 — T / 
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• If v\ > V 2 , then 


r {X2-V2 X2~V2 + VJ\ (X2-V2 + W X2~V2\ r W 

Ia,S\ -,- = Ia,s[ -,- + d — 

\ T T / \ T T / T 


Notice that 




1 — r / VI — r/ V 1 — r 
Since Ia is convex and v\ > u 2 , we have 


1 — t J 1 — T 


ia( v 4^] - - // 1,1 


1 — T 


1 — r 


1 — r 


1 — r 


Thus, 


i s \ v 4^± ) + OA-1 > is(A^-) +uf^i 9 


1 — r/ VI —r/ Vl — r/ V 1 - r / 1 —r 

Since and are convex and X 2 — V 2 + w + Vi — Q — we have 

5 ( 9 ) > ( 5 w + /a,s(xi,x 2 ) . 

We conclude that J(xi,X 2 ,w) = 5w + Ia,s(x i,x 2 ). 


□ 


In other words, the above theorem states that if the input variables satisfy large deviations 
principles with rate functions I a and Is with va the d'-tilted measure of i/g, then the output 
variables satisfy large deviations principles with the same rate functions. 


References 

[1] Burke, P., The output of a queueing system. Operations Research , 4:699-704, 1956. 

[2] Chang, C., Stability, queue length and delay of deterministic and stochastic queueing 
networks. IEEE Transactions on Automatic Control, 39:913-931, 1994. 

[3] Courcoubetis, C., Kesidis, G., Ridder, A., Walrand, J. and Weber, R., Admission control 
and routing in ATM networks using inference from measured buffer occupancy. IEEE 
Transations on Communications, 43:1778-1784, 1995. 

[4] Dernbo, A. and Zajic, T., Large deviations : from empirical mean and measure to partial 
processes. Stochastic Processes and their Applications, 57:197-224, 1995. 

[5] Dernbo, A. and Zeitouni, O., Large Deviations Techniques and Applications. Jones and 
Barlett, Boston, 1993. 

[6] Draief, M., Mairesse, J. and O’Connell, N., Joint Burke’s theorem and RSK representation 
for a queue and a store. In C. Banderier and C. Krattenthaler, editors, Discrete Random 
Walks, DRW'03, volume AC, pages 69-82, Institut Henri Poincare, Paris, 2003. 

[7] Draief, M., Mairesse, J. and O’Connell, N., Queues, Stores, and Tableaux. To appear in 
Journal of Applied Probability 42(4), 2005. 

[8] Duffield, N. and O’Connell, N., Large deviations and overflow probabilities for the gen¬ 
eral single server queue, with applications. Mathematical Proceedings of the Cambridge 
Philosophical Society, 118(1), 1995. 


19 



[9] Gallager, R. and Prabhakar, B., Entropy and the timing capacity of discrete queues. IEEE 
Transactions on Information Theory, 49(2):357-370, 2003. 

[10] Ganesh, A. and O’Connell, N., A large deviation principle with queueing applications. 
Stochastics and Stochastic Reports, 73(l-2):25-35, 2002. 

[11] Ganesh, A., O’Connell, N. and Prabhakar, B., Invariant rate functions for discrete-time 
queues. Annals of Applied Probability, 13(2):446-474, 2003. 

[12] Ganesh, A., O’Connell, N. and Wischik, D., Big Queues. Lecture Notes in Mathematics. 
Springer-Verlag, Berlin, 2004. 

[13] Glynn, P. and Whitt, W., Logarithmic asymptotics for steady-state tail probabilities in a 
single-server queue. Journal of Applied Probability, 31A:131-156, 1994. 

[14] Kelly, F., Notes on effective bandwidth. In F. Kelly and I. Ziedens, editors, Stochastic 
Networks, pages 141-168. Oxford University Press, 1996. 

[15] Kelly, F., Effective bandwidths at multi-class queues. Queueing Systems Theory and 
Application, 9:5-15, 1991. 

[16] Loynes, R., The stability of a queue with non-independent interarrival and service times. 
Mathematical Proceedings of the Cambridge Philosophical Society, 58:497-520, 1962. 


20 



